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ABSTRACT 

We present a numerical study of the dynamics of magnetized, relativistic, non-self-gravitating, 
axisymmetric tori orbiting in the background spacetimes of Schwarzschild and Kerr black 
holes. The initial models have a constant specific angular momentum and are built with a non- 
zero toroidal magnetic field component, for which equilibrium configurations have recently 
been obtained. In this work we extend our previous investigations which dealt with purely hy- 
drodynamical thick discs, and study the dynamics of magnetized tori subject to perturbations 
which, for the values of the magnetic field strength considered here, trigger quasi-periodic 
oscillations lasting for tens of orbital periods. Overall, we have found that the dynamics of 
the magnetized tori analyzed is very similar to that found in the corresponding unmagnetized 
models. The spectral distribution of the eigenfrequencies of oscillation shows the presence of 
a fundamental p mode and of a series of overtones in a harmonic ratio 2:3: . . . . These sim- 
ulations, therefore, extend the validity of the model of iRezzoUa et al.l ( l2003a[) for explaining 
the high-frequency QPOs observed in the spectra of LMXBs containing a black-hole candi- 
date also to the case of magnetized discs with purely toroidal magnetic field distribution. If 
sufficiently compact and massive, these oscillations can also lead to the emission of intense 
gravitational radiation which is potentially detectable for sources within the Galaxy. 

Key words: accretion discs - general relativity - hydrodynamics - oscillations - gravitational 
waves 



1 INTRODUCTION 



In a series of recent pap ers jZanotti et alj l2003l : iRezzolla et al.l 
l2003bl :l Zanotti et al.ll2005l) it has been shown that upon the intro- 
duction of perturbations, stable relativistic tori (or thick accretion 
discs) manifest a long-term oscillatory behaviour lasting for tens 
of orbital periods. When the average disc density is close to nu- 
clear matter density, the associated changes in the mass-quadrupole 
moment make these objects promising sources of high-frequency, 
detectable gravitational radiation for ground-based interferometers 
and advanced resonant bar detectors, particularly for Galactic sys- 
tems. This situation applies to astrophysical thick accretion discs 
formed following binary neutron star coalescence or the gravita- 
tional core collapse of a sufficiently massive star. If the discs are 
instead composed of low-density material stripped from the sec- 
ondary star in low-mass X-ray binaries (LMXBs), their oscilla- 
tions could help explaining the high-frequency quasi-periodic os- 
cillations (QPOs) observed in the spectra of X-ray binaries. In- 
deed, such QPOs can be explained in terms of p-mode oscilla- 
tions of_a_small2si^e_torus orbiting around a stellar-mass black 
hole jRezzoUa et alj2003ah . 

The studies reported in the papers mentioned above have con- 



sidered both Schwarzschild and Kerr black holes as well as constant 
and nonconstant (power-law) distributions of the specific angular 
momentum of the discs. However, they have so far been limited 
to purely hydrodynamical matter models, neglecting a fundamen- 
tal aspect of such objects, namely the existence of magnetic fields. 
There is general agreement that magnetic fields are bound to play 
an important role in the dynamics of accretion discs orbiting around 
black holes. They can be the sou rce of viscous processes w ithin the 
disc through MHD-turbulence dShakura & Sunvaev|[l973l) , as con- 
firmed by t he presence of the so called magnetorotational instabil- 
ity (MRI) tealbusi 12003')) that regulates the accretion process by 
transferring angular momentum outwards. In addition, the forma- 
tion and collimation of the strong relativistic outflows or jets rou- 
tinely observed in a variety of scales in astrophysics (from micro- 
quasars to radio-galaxies and quasars) is closely linked to the pres- 
ence of magnetic fields. 

General relativistic magnetohydrodynamic (GRMHD here- 
after) numerical simulations provide the best approach for the 
investigation of the dynamics of relativistic, magnetized accre- 
tion discs under generic nonlinear conditions. In recent years 
there have been important breakthroughs and a sustained level 
of activity in the modelling of such systems, as formula- 
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tions of the GRMHD equations in forms suitable for numeri- 
cal work have become available. This has been naturally fol- 
lowed by their implementation in state-of-t he-art numerical codes 
develo p ed by a number of groups {e.g., lOe Villier s & HawlevI 
(l2003h- iGamm ie. McKinnev & Toth (2003); Komissaro TfeOOSh: 
iDuez et aljlToOSi): iShibata & Sekiguchii t20 05): Anni nos et alj 



ll2005h: iFragild 12005 ): 

iMizuno et al. 1 20061) : Giacomazzo & Rezzolla 



Anton et al. 



l200a)~ lMcKinne\l ( l2006h : 
(2007)) many of 
which have been applied to the investigation of issues such as 
the M RI in accreti o n disc s and jet formation. Moreover, very re- 
cently iKomissarovl ( |2006|) has derived an analytic solution for an 
axisymmetric, stationary torus with constant distribution of specific 
angular momentum and a toroidal magnetic field configuration that 
generalizes to the relativis t ic reg ime a previous Newtonian solu- 
tion found bv lOkada et aP ( Il989l) . Such equilibrium solution can 
be used not only as a test for GRMHD codes in strong gravity, but 
also as initial data for numerical studies of the dynamics of magne- 
tized tori when subject to small perturbations. The latter is, indeed, 
the main purpose of the present paper. 

In this way we aim at investigating if and how the dynam- 
ics of such objects changes when the influence of a toroidal mag- 
netic field is taken into account. We discuss the implications of our 
findings on the QPOs observed in LMXBs with a bl ack hole candi- 
date, as sessing the validity of the model proposed bv lRezzoUa et al.l 
( l2003ah in a more general context. 

The paper is organized as follows: in Section 2 we briefly re- 
view the equilibrium solution found by Komissarov (2006) for a 
stationary torus with a toroidal magnetic field orbiting around a 
black hole. The mathematical framework we use for the formu- 
lation of the GRMHD equations and for their implementation in 
our numerical code is discussed in Section 3, while in Section 4 
we describe the approach we follow for the numerical solution of 
the GRMHD equations. Section 5 is devoted to the discussion of 
the initial models considered, with the results being presented in 
Section 6. Finally, Section 7 summarizes the paper and our main 
findings. We adopt a geometrized system of units extended to elec- 
tromagnetic quantities by setting G = c = eo = 1. where eo is 
the vacuum permittivity. Greek indices run from to 3 and Latin 
indices from 1 to 3. 



2 STATIONARY FLUID CONFIGURATIONS WITH A 
TOROIDAL MAGNETIC FIELD 

The initial configurations we consider can be considered as 
the MHD extensions to of the stationary hydrodynamical solu- 
tions of thick discs o r biting around a black ho le described by 
iKozlowski et alj ( ll978l) . lAbramowicz et al] ( Il978l) and are built us - 
ing the analytic solution suggested recently by Komissarovl J2006h . 

The basic equations that are solved to construct such initial 
models are the continuity equation V ^{pu^) = for the rest-mass 
density p, the conservation of energy-momentum V^jT*^" = 0, and 
Maxwell's equation V F'^'^) — 0, where the operator is the 
covariant derivative with respect to the spacetime four-metric and 
*F'^'^ is the dual of the Faraday tensor defined as 



F — u b 



u b 



(1) 



In this expression is the fluid four-velocity and b^ is the mag- 
netic field measured by an observer comoving with the fluid. As 
usual in ideal relativistic MHD (i.e., for a plasma having infinite 



conductivity), the stress-energy tensor T^" is expressed as 

T"-^ = {ph + b^)u^u-' +fp+^\ _ b^b-, 



(2) 



where g''" are the metric coefficients, p is the (thermal) pressure, h 
the specific enthalpy, and 6^ = b'^b^. 

The equilibrium equations are then solved to build station- 
ary and axisymmetric fluid configurations with a toroidal magnetic 
field distribution in the tori and a constant distribution of the spe- 
cific angular momentum in the equatorial plane. T he main differ- 
ence of our solution with that o f IKomissarovl ilOOdt) is that we em- 
ploy a polytropic equation of state (EOS) of the form p = Kp^ for 
the fluid, where k is the polytropic constant and F is the adiabatic 
index. Such an EOS has a well-defined phys ical meaning and dif- 
fers from the one used bv lKomissarovl j2006l) , p = Kui'', where uj 
is the fluid enthalpy, and K and q are constants. 

By imposing the condition of axisymmetry and stationarity in 
a spherical coordinate system (i.e., = dt = 0), the hydrostatic 
equilibrium conditions in the r and 6 directions are given by 



Vi \n{ut) - + + -— — 

1 — i2£ vj 2Lw 



= 0, 



(3) 



with i ^r,6 and C{r, 6) = gt4,gtci, 
appearing in ([3} is defined as 



g<t,<pgtt. The angular velocity 



ll'' 



the specific angular momentum is given by 

Ut 

and the components of the magnetic field are 



(4) 



(5) 



b*" = ^2p™/(g^^ + 243t^ + A)ff«), (6) 

6' = lob''' . (7) 

Following' Komissarovl ( |200( ^. we consider the following EOS 
for the magnetic pressure Pm = MC'^~^w'^ , where M and q are 
constants, and which essentially amotmts to confining the magnetic 
field to the interior of the torus. Using this relation, we can integrate 
eq. (|3]l, which in the case of constant specific angular momentum 
yields 



W-Win + ln 1 



VK 



+ - 



-M {Cw 



^ 0, (8) 



where the potential W is defined as = In | . Note that in gen- 
eral there will be two radial locations at which io equals the Ke- 
plerian specific angular momentum. The innermost of these radii 
represents the location of the "cusp" of the torus, while the outer- 
most the "centre". When a magnetic field is present, the position of 
the centre does not necessarily correspond with that of the pressure 
maximum, as in the purely hydrodynamical case. 

In order to solve Eq. ([3} a number of parameters are needed to 
define the initial model, namely F, q, £o. Win and the ratio of the 
magnetic-to-gas pressure at the centre of the torus, Pc = {pm/p)c. 
Thus, using the definition of /9c, we obtain the rest-mass density at 
the centre of the torus from the following expression: 



Wc - Win + In (^H 
/9cF 



F - 1 



kpc ^ 



(r- 



F/(F-1) 



(9) 



Finally, the equilibrium equation ([3} can be solved to obtain 
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the distributio n of all relevant mag netohydrodynamic quantities in- 
side the torus iKomissarovl ( l2006h ). 



3 GENERAL RELATIVISTIC MHD EQUATIONS 

As mentioned in the Introduction, there has been intense work 
in recent years on formulations of the GRMHD equations suit- 
able for numerical app roaches ( Gammie, McKinney & Toth 2003 ; 
De Villiers & Hawley 2003; Komissarov 2005; Duez et al. 2005; 
Shibata & Sekiguch i 2005; Ann inos et a l. 2005; Anton et al. 2006. ; 
Giacomazzo & Rezzolla 2007). We here follow the approach laid 



Anton et^ |2006) and adopt the 3 + 1 formulation of gen- 



out in 

era! relativity in which the 4-dimensional spacetime is foliated into 
a set of non-intersecting spacelike hypersurfaces. The 3 + 1 line 
element of the metric then reads 



ds 



'{a ~ f3il3')dx°dx° + 2/3idx'dx° + jijdx'dx^ , (10) 



where 7ij is the 3-metric induced on each spacelike slice, and a 
and /3' are the so-called lapse function and shift vector, respectively. 
Under the ideal MHD condition. Maxwell's equations 
_ g j.gjjjjj.g f[jg divergence-free condition for the mag- 
netic field 



0, 



(11) 



together with the induction equation for the evolution of the mag- 
netic field 

d 



1 d . _ 1 

-X^B ) = 



V^dx- 



■{^[av'B' - aS'B']}, (12) 



where 7 = det(7ij) and ff' = v' — /ct, with and being 
respectively the spatial components of the velocity and of the mag- 
netic field, as measured by the Eulerian observer associated to the 

3 + 1 splitting. 

Following I Anton et alj ( 1200^ , the conservation equations for 
the energy-momentum tensor given by Eq. ^ together with the 
continuity equation and the induction equation for the magnetic 
field can be written as a first-order, flux-conservative, hyperbolic 
system. The state vector and the vector of fluxes of the fundamen- 
tal GRMHD system of equations read 



1 / 9V7F° ^ d 

/—g \ dx'^ dx 



S, 



(13) 



where g = det((7^^) = 0^7. The state vector F" is given by 



D 

r 



with the definitions 

D = pW, 

Sj = p{h + b'^/p)W^Vj - ab°bj , 
T = p{h + h''lp)W^ -(p + b^/2)~o?(b"f 



(14) 



(15) 
(16) 
D, (17) 



and where W is the Lorentz factor of the fluid. The "fluxes" F* in 
eqs. J13t have instead explicit components given by 



F' 



Sjv' + (p + b'^/2)S] - bjB'/W 
ri,' + {p + b'^/2)v' - ab°B^/W 



(18) 



while the "source" terms S are 


T^'- {dg.,/dx'' - 
a {T^'^dlna/dx" 



rpjJ.l'-pO 



(19) 



where O'' = (0,0,0)^, and F'^g are the Christoffel symbols for 
either a Schwarzschild or Kerr black-hole spacetime. Note that the 
following fundamental relations hold between the four components 
of the magnetic field in the comoving frame, b'^, and the three vec- 
tor components B^ measured by the Eulerian observer 



,0 WB'v, 
b — 



b' = 



B' + ab\' 
W 



(20) 



(21) 



(22) 



Finally, the modulus of the magnetic field can be written as 

2_ B^+aW 

where = B'B,. 

Casting the system of evolution equations in flux- 
conservative, hyperbolic form allows us to take advantage of 
high-resolution shock-capturing (HRSC) methods for their numer- 
ical solution. The hyperbolic structure of those equations and the 
associated spectral decomposition of the flux-vector Jacobians, 
n eeded for their num erical solution with Riemann solvers, is given 
in lAnt6netal] ( l2006l) . 



4 NUMERICAL APPROACH 

The numerical code used for the simulations report ed in this pa- 
per is an ext ended version of the code presented in I Zanotti et al.l 
(20031 l2005h to account for solution of the GRMH D equations. 
The ac curacy of the code has been recently assessed in lAnton et al.l 
1 I2OO6I) , with a number of tests including magnetized shock tubes 
and accretion onto Schwarzschild and Kerr black holes. The system 
of GRMHD equations J13t is solved using a conservative HRSC 
scheme based on the HLLE solver, except for the induction equa- 
tion for which we use th e co nstraint transport method designed by 
lEvans & Hawley| l ll988l) and lRvu et alj ( ll998l) . Second-order accu- 
racy in both space and time is achieved by adopting a piecewise- 
linear cell reconstruction procedure and a second-order, conserva- 
tive Runge-Kutta scheme, respectively. 

The code makes use of polar spherical coordinates in the two 
spatial dimensions (r, 6) and the computational grid consists of 
Nr X Ne zones in the radial and angular directions, respectively. 
The innermost zone of the radial grid is placed at rmin = ''horizon + 
0.1, and the outer boundary in the radial direction is at a distance 
about 30% larger than the outer radius of the torus, rout - The radial 
grid has typically Nr — 300 and is built by joining smoothly a first 
patch which extends from rmin to the outer radius of the torus and 
is logarithmically spaced (with a maximum radial resolution at the 
innermost grid zone, Ar/M = 1 x 10~^, where M is the mass of 
the black hole) and a second patch with a uniform grid and which 
extends up to rmax. On the other hand, the angular grid consists of 
Ng — 100 equally spaced zones and covers the domain from to 

TV. 

As in the hydrodynamical code, a low density atmosphere is 
introduced in those parts of the computational domain not occu- 
pied by the torus. This is set to follow the spherically-symmetric 
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Table 1. From left to right the columns report the name of the model, the spin of the black hole, a, the specific angular momentum, Iq, the polytropic constant, 
K, the inner and outer radius of the torus, and r'out, the orbital period at the point of maximum rest-mass density, torb- the maximum rest-mass density, 
Pmax, the magnetic-to-gas pressure at the maximum of the rest-mass density, /3c, and the maximum magnetic field, -Bmax- For all models the torus-to-hole 
mass ratio Mt/M is 0.1 (M = 2.5M0) and the adiabatic exponent of the EOS is 4/3. 



Model 


a 




K (Cgs) 




Tout 


tovh (ms) 


Pmax (cgs) 


ft 


Bmax (G) 


51 


0.0 


3.80 


9.33x10^^ 


4.57 


15.88 


1.86 


1.25x10" 


0.00 


0.0 


52 


0.0 


3.80 


9.21 xlO^^ 


4.57 


15.88 


1.86 


1.26x10" 


0.01 


2.50x10" 


53 


0.0 


3.80 


9.10x10^3 


4.57 


15.88 


1.86 


1.27x10" 


0.02 


3.52x10" 


54 


0.0 


3.80 


8.90x10^3 


4.57 


15.88 


1.86 


1.28x10" 


0.04 


4.94x10" 


55 


0.0 


3.80 


8.40x10^3 


4.57 


15.88 


1.86 


1.29x10" 


0.10 


7.58x10" 


56 


0.0 


3.80 


7.60x10^3 


4.57 


15.88 


1.86 


1.34x10" 


0.20 


1.04x10" 


57 


0.0 


3.80 


6.00x10^3 


4.57 


15.88 


1.86 


1.39x10" 


0.50 


1.50x10" 


58 


0.0 


3.80 


4.49x10^3 


4.57 


15.88 


1.86 


1.40x10" 


1.00 


1.85x10" 


Kl 


0.5 


3.30 


2.20x10^'' 


3.16 


15.65 


1.22 


1.44x10" 


0.01 


4.29x10" 


K2 


0.7 


3.00 


2.25 xlO^'' 


2.57 


12.07 


0.88 


2.74x10" 


0.01 


6.69x10" 


K3 


0.9 


2.60 


7.80x10^'' 


1.77 


19.25 


0.56 


1.87x10" 


0.01 


1.15x10" 



accreting solution described bv lMichell ( [T972h in the case that the 
background metric is that of a Schwarzschild black hole and a mod- 
ified solution, whic h accounts for the rotation of the black hole 
( IZanottietal] ( l2005h ), when we consider the Kerr background met- 
ric. Since this atmosphere is evolved as the rest of the fluid and is 
essentially stationary but close to the torus, it is sufficient to ensure 
that its dynamics does not affect that of the torus. This is the case if 
the maximum density of the atmosphere is 5-6 orders of magnitude 
smaller than the central density of the torus. Note that since we limit 
our analysis to isoentropic evolutions of isoentropic initial models, 
the energy equation needs not to be solved. Finall y, the boundary 
condit ions adopted are the same as those used by iFont & Daignd 
( l2002h . 



5 INITIAL MODELS 

The initial models consist of a number of magnetized relativistic 
tori which fill their outermost closed equipotential surface, so that 
their inner radii coincide with the position of the cusp, rin = ^cusp- 
In practice, we determine the positions of the cusp and of the max- 
imum rest-mass density in the torus by imposing that the specific 
angular momentum at these two points coincides with the Keple- 
rian value. Clearly, different values of specific angular momentum 
will produce tori with different positions of the cusp and of the 
maximum rest-mass density. In a purely hydrodynamical context, 
the effect on the dynamics of the tori of the distribution of specific 
angular momentum, bei ng either cons t ant or s atisfy ing a power-law 
with r, was studied by IZanotti et al.l ( |2003| , l2005h . In this paper, 
however, we consider only magnetized tori with constant specific 
angular momentum as we want to first focus on the influence a 
magnetic field has on the dynamics, both in a Schwarzschild and in 
a Kerr background metric. In this way we can conveniently exploit 
the analytic solution reviewed in Sect. 2 and which cannot be ex- 
tended simply to include the case of nonconstant specific angular 
momentum distributions. 

Once the specific angular momentum is fixed, the inner edge 
of the torus rin is determined by the potential gap at such inner 
edge, A Win = Win — Wcusp which, in the case of constant specific 
angular momentum distributions, is defined as 

AWin = ln[(-ltt)in] - ln[(-t(t)cuBp], (23) 



with AWin — corresponding to a torus filling its outermost 
equipotential surface. 

All of the models are built with an adiabatic index F = 4/3 to 
mimic a degenerate relativistic electron gas, and the polytropic con- 
stant K is fixed such that the torus-to-black hole mass ratio, Mt/M, 
is roughly 0.1. Since the mass of the torus is at most 10% of that 
of the black hole, we can neglect the self-gravity of the torus and 
study the dynamics of such objects in a fixed background space- 
time (test-fluid approximation). Moreover, the disc-to-hole mass 
ratio adopted here is in agreement with the one obtained in sim- 
ul ations of unequal m ass binary neutron star mergers performed 
bv lShibata et"dT ( l2003l) and lShibata & Sekiguchil ( |2005|) . 

Overall, we have investigated a number of different models for 
tori orbiting either nonrotating or rotating black holes. In the case of 
Schwarzschild black holes, the main difference among the models 
is the strength of the toroidal magnetic field, which is parametrized 
by the ratio of the magnetic-to-gas pressure at the centre of the 
disc, /3c = bV(2p)- In the case of Kerr black holes, on the other 
hand, we report results for tori orbiting around black holes with 
spins a — 0.5, 0.7 and 0.9, while keeping constant the magnetic- 
to-gas pressure ratio at /3c = 0.01. A summary of all the models 
considered is given in Table[T] 

The set of models chosen here will serve a double purpose. Be- 
ing tori with large average densities, they can provide accurate esti- 
mates for the gravitational-wave emission triggered by the oscilla- 
tions. On the other hand, since the ratio among the eigenfrequencies 
is the astrophysically most relevant quantity and this does not de- 
pend on the density, this set of models is also useful for analysing 
the oscillation properties of the accretion discs in LMXBs. It is also 
important to note that for tori with /3c > 1 the initial solution de- 
grades over time as a significant mass is accreted in these cases, 
with an accretion rate that increases with the strength of the mag- 
netic field. The dependence of the stability of thick discs with the 
strength of the toro idal magnetic fi eld will be the subject of an ac- 
companying paper teezzolla et all ( 120071) ). 

The maximum strength of the magnetic field at the centre, de- 
termined by the parameter /3c, can be calculated through Eq. 
which also reflects the dependence of the toroidal magnetic field 
component on the background metric. The initial models consid- 
ered are such that /3c takes values between and 1, as shown 
in Table [T] This also fixes the overall strength of the magnetic 
field, whose maximum values are reported in the tenth column 
of the same table. The values of the magnetic field strength at 
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Figure 1. Isocontours of the logarithm of the rest-mass density for the unperturbed model 52. The left panel shows the configuration at the initial time and the 
light panel the con'esponding distribution after 100 orbital timescales. The equilibrium solution is preserved to high accuracy. 



the centre for the case of tori around a Schwarzschild black hole 
range from the magnetized model S2 with 2.50x10^^ Gauss to 
model S8 with 1.85x 10^® Gauss. These values are in good agree- 
ment with the typical values expected to be present in the as- 
trophysical scenarios that could for m a relativistic thick torus , 
such as the ma gnetized core collapse (Cerda-Duran & Fontll2006l ; 
lObergaulinger et al. 2006; Shibata et al. 2006) and values consid- 
ered for the collapse of magnetized hypermassive neutron stars by 
Duez et al. (2006). 

In order to trigger the oscillations, we perturb the models re- 
ported in Table [T] by adding a small radial velocity (we recall that 
in equilibrium all velocity compo nents but the azimuthal on e are 
zero). As in our previous work jZanotti et al.l |2003| , 1 20051) , this 
perturbation is parametrized in terms of a dimensionless coeffi- 
cient 7? of the sphe rically symmetric accretion flow on to a black 
hole ( lMiche]|[l973) , i.e., = r;(w,.) Michel- In all the simulations 
reported we choose = 0.1, but the results are not sensitive to 
this cho ice as long as t he osc illations are in a linear regime {i.e., for 
?7 < 0.2 lZanotti et"ZH2003h ). 

Finally, it is worth commenting on the choice made for the 
magnetic field distribution. On the one hand this choice is mo- 
tivated by the mere convenience of having an analytic equilib- 
rium solution upon which a perturbation can be introduced. On 
the other hand, there exists an additional motivation which is more 
astrophysically motivated. As it has been shown in recent simu- 
l ations of m a gnetized co r e collapse (C erda-Duran & Font 2006 ; 
[bbergauUnger et alj|2003 : IShibata et alJl2006D the magnetic field 
distribution in the nascent, magnetized, proto-neutron stars has a 
dominant toroidal component, quite irrespective of the initial con- 
figuration. Since gravitational core collapse is one of the processes 
through which thick accretion discs may form, the toroidal initial 
configuration of our simulations is well justified. This choice, how- 
ever, also has an important consequence. Because of the absence of 
an initial poloidal magnetic field, in fact, the magneto-rotational in- 
stability (MRI), which could change even significantly the dynam- 



ics of our tori, cannot develop in our simulations. Indeed, iFragild 
(2005) has investigated the oscillation of an accretion torus hav- 
ing an initial poloidal magnetic field component. Although pre- 
liminary, his results suggest that the development of the MRI an d 
of the Papaloizou-Pringle instability dPapaloizou &Pringlel ( IT984l) ) 
may damp significantly the oscillation modes of accretion tori with 
poloidal magnetic fields. We will address this question in a future 
work. 



6 RESULTS 

6.1 Oscillation properties 

6. 1. 1 Dynamics of magnetized tori 

We have first investigated equilibrium configurations of magne- 
tized tori by performing numerical evolutions of unperturbed tori 
(r) — 0) and by checking the stationarity of the solution over a 
timescale which is a couple of orders of magnitude larger that the 
dynamical one. As representative example, we show in Fig. [T] the 
isocontours of the logarithm of the rest-mass density of model S2 
as computed at the initial time t = (left panel) and at the time 
when the simulation was stopped (right panel). This corresponds 
tot = 100 torh, where torh is the Keplerian orbital time for a par- 
ticle in a circular orbit at the centre of the torus. Aside from the 
minute accretion of matter from the cusp towards the black hole 
(see below), the final snapshot of the rest-mass distribution clearly 
shows the stationarity of the equilibrium initial solution. More pre- 
cisely, the central rest-mass density, after a short initial transient 
phase, settles down to a stationary value which differs after 100 
orbital timescales only of 2% from the initial one. This provides 
a strong evidence of the ability of the code to keep the torus in 
equilibrium for evolutions much longer than the characteristic dy- 
namical timescales of these objects. 

On the left panel of Fig.|2]we show instead the evolution over 
lOOforb of the central rest-mass density of the least magnetized 
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model 52, when a perturbation with parameter -q = 0.1 is added 
to the equilibrium modeQ. It is interesting to note that despite the 
presence of a rather strong toroidal magnetic field, the persistent 
oscillatory behaviour found in these simulatio ns is very similar to 
the o ne found in purely hydrodynamical tori jZanotti et aLlbOOl 
1200 5). Note also that the small secular decrease in the oscilla- 
tion amplitude is not to be related to numerical or physical dissi- 
pation, since the code is essentially inviscid and the EOS used is 
isoentropic. Rather, we believe it to be the result of the small but 
nonzero mass spilled through the cusp at each oscillation (see also 
discussion below). Furthermore, on a smaller timescale than the 
one shown in Fig. [2] the oscillations show a remarkable harmonic 
behaviour and this is highlighted in th e small insets i n Fig.|2l This 
is in stark contrast with the results o f lFragild j2005l) . which were 
obtained with comparable numerical resolutions, but with an ini- 
tial poloidal magnetic field configuration. In that case, in fact, the 
oscillations were rapidly damped in only a few orbital periods. 

Results from a representative model with a higher magnetic 
field are shown in the right panel of Fig. |2l which again reports 
the evolution of the normalized central rest-mass density for model 
58. Note that despite this model has a magnetic-to-gas pressure 
ratio at the centre Pc = 1, and hence a central magnetic field of 
~ 2 X 10^^ G, its overall the dynamics is very similar to that of 
model S2. Also in this case, in fact, the oscillations are persistent 
during the entire evolution (100 torb) and show almost no damping. 
However, the amplitude does show variations over time and, most 
importantly, it no longer maintains a symmetric behaviour between 
maxima and minima, as a result, we believe, of the increased mass 
accretion through the cusp. We recall, in fact, that all the initial 
models considered in our sample correspond to marginally stable 
tori, i.e., tori filling entirely their outermost closed equipotential 
surface. Any perturbation, however small, will induce some mat- 
ter to leave the equipotential surface through the cusp, leading to 



^ We have here chosen to show the evolution of the rest-mass density as 
this is has a simple physical interpretation, but all of the MHD variables 
exhibit the same harmonic behaviour. 



the accretion of mass and angular momentum onto the black hole. 
Evidence in favour of this is shown in Fig. [3] which reports the ac- 
cretion mass-flux for model S2 (upper panel) and 58 (lower panel). 
While both reflect the oscillations in the dynamics, they also have 
different mean values, with the one relative to model 58 being al- 
most an order of magnitude larger. Note also the correlation be- 
tween the fluctuations in the mass-accretion rate and the changes in 
the oscillation amplitudes shown in Fig.|2] In particular, the sudden 
change in the mass-flux of model 58 at f ~ 35 torb and which cor- 
responds to a change in the amplitude modulation in the right panel 
of Fig. [2] 

Although the accretion-rates are well above the Eddington 
limit (which is ~ 10"^'' Mq/s for a 2.5 Mq black hole), the 
amounts of mass accreted by the black hole at t = 100 torb is only 
1.3% and 3.3% of the initial mass for models 5*2 and 5*8, respec- 
tively. Similarly, the total amount of angular momentum accreted 
at the end of the simulation would introduce a change in the black- 
hole's spin of less than 1% for both models S2 and S8. Overall, 
therefore, these changes in the mass and spin of the black holes 
are extremely small and thus justify the use of a fixed background 
spacetime. Finally, in Fig.|4]we show the evolution of the normal- 
ized central rest-mass density for model K2, which corresponds to 
a torus orbiting around a Kerr black hole with spin a = 0.7. Again, 
a perturbation with parameter rj — . 1 was added to the equilibrium 
model so as to investigate the oscillatory behaviour of the torus 
around its equilibrium position. As in the purely hydrodynamical 
case, the qualitative behaviour in models around a Kerr black hole 
is very similar to that found for models around a Schwarzschild 
black hole, and the dynamics shows, also in this case, a negligible 
damping of the oscillations after the initial transient. 



6.1.2 Power Spectra 

An important feature of axisymmetric p-mode oscillations of ac- 
cretion tori is that the lowest-order eigenfrequencies appear in the 
harmonic sequence 2 : 3. This feature was first discovered in 
the purely hydrodynamical numerical simulations of IZanotti et al.l 
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Figure 3. Time evolution of the mass accretion rate for models 52 and S8. 

( I2OO3I) , subsequently confirmed through a perturbative analysis 
in a Schwarzschild spacetime by Rezzollaetal. (2003b), and 
later extended to a Kerr spacetime and to more g eneral d i stribu - 
tion s of the specific ang ular momentum by iZanotti et al. I l l2005h 
and lMontero et alj ( |2004|) . Overall, it was found that the 2 : 3 har- 
monic sequence was present with a variance of ~ 10% for tori with 
a constant distribution of specific angular momentum and with a 
variance of ~ 20% for tori with a power-law distribution of spe- 
cific angular momentum. Since the 2 : 3 harmonic sequence is the 
result of global modes of oscillation, it depends on a number of 
different elements that contribute to small deviations from an exact 
relation among integers. The latter, in fact, should be expected only 
for a perfect one-dimensional cavity, trapping the p modes without 
losses. In practice, however, factors such as the vertical size of the 
tori, the black hole spin, the distribution of specific angular momen- 
tum, the EOS considered, and the presence of a small but nonzero 
mass-loss, can all influence this departure. 

While the understanding of the properties of these modes 
of o scillation has grown considerably over the last few years 
(see iMontero et al.l (2004) for a list of references), and an ex- 
haustive analysis ha s been made in the case of relativistic slen- 
der tori qBlaes et al. I l l2006h ), it was not obvious whether such a 
harmonic sequence would still be present in the case of magne- 
tized discs with toroidal magnetic fields. To address this ques- 
tion, we have performed a Fourier analysis of the time evolu- 
tion of some representative variables and obtained quantitative in- 
formation on the quasi-periodic behavior of the tori. In particu- 
lar, for all of the models considered, we have Fourier-transformed 
the evolution of the L2 norm of the rest-mass density, defined as 
l|p||2 = Si^Ti X]jlfi(Pij)^ studied the properties of result- 
ing power spectra. These, we recall, show distinctive peaks at the 
frequencies that can be identified with the quasi-normal modes of 
oscillation of the disc0. Clearly, the accuracy in calculating these 
eigenfrequencies depends linearly on the length of the timeseries 

^ Note that because of the underlining axisymmetry of our calculations, 
we cannot compute the effect of transverse hydromagnetic waves, such as 
Alfven waves, propagating along the toroidal magnetic field lines. 



Figure 4. Time evolution of the central rest-mass density normalized to its 
initial value for the model K2 



and is of 0.01 kHz for the evolutions carried out here and that ex- 
tend for lOOtorb ~ 100 ms. 

In Fig.|5]we present the power spectra (PSD) obtained from 
the L2 norm of the rest-mass density for model S2 (left panel) 
and 58 (right panel); in both panels the solid lines refer to the 
magnetized tori, while the dashed ones to the unmagnetized coun- 
terpart SI, which is shown for reference. A rapid look at the 
panels in Fig. |5] reveals that the overall dynamics of magnetized 
tori shows features whi c h are surprisingly similar to those found 
by IZanotti et al. I l l200l I2OO5I) for unmagnetized accretion tori. 
Namely, the spectra have a fundamental mode / (which is the mag- 
netic equivalent of the p mode discussed in IZanotti et alj ( l2003l 
l2005h ) and a series of overtones, for which, in particular, the first 
overtone oi can usually be identified clearly. Interestingly, also 
these spectra show the 2 : 3 harmonic relation between the fre- 
quencies of the fundamental mode and its first overtones. Such a 
feature remains therefore unmodified and an important signature of 
the oscillations properties of magnetized tori with a toroidal mag- 
netic field. 

It is also worth noting that in the case of mildly magnetized 
tori, such as model S2, the similarity in the PSD is rather striking 
and the two spectra differ only in the relative amplitude between the 
eigenfrequencies 01,02,... and the modes which are the result of 
nonlinear coupling {e.g., 2/ — oi , 2/, . . .). On the other hand, in the 
case of more highly magnetized tori, such as model 58, the mag- 
netic field strength is sufficiently large to produce variations in the 
eigenfrequencies, which are all shifted to higher frequencies, with 
deviations, however, which become larger for higher overtones. 
While not totally unexpected [a magnetic field is know n to in- 
crease the eigenfrequencies of magnetized stars dNasiri & Soboutil 
([198^))], these represent the first calculations of the eigenfrequen- 
cies of relativistic magnetized discs and, as such, anticipate analo- 
gous perturbative studies. 

As a way to quantify the differential shift of the eigenfrequen- 
cies to larger values, we report in Table [2] the frequencies of the 
fundamental mode, of the first overtone, and their ratio for all of the 
models considered. Another analogy worth noticing in the spectra 
presented in Fig. [5] is the presence of nonlinear couplings among 
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Figure 5. Solid lines: Power spectram of the L-2 norm of the rest-mass density for models 5*2 (left panel) and SS (right panel). The dashed lines show the 
corresponding PSD for the unmagnetized versions of both models. The units in the vertical axis are arbitrary and the PSDs were obtained using a Manning 
filter 



Table 2. From left to right, the columns report the name of the model, the 
frequency of the fundamental mode, the frequency of the first overtone, their 
ratio, and the magnetic-to-gas pressure ratio at the centre of the torus. 



Model 


/(Hz) 


oi (Hz) 


oi/f 


/3c 


SI 


224 


332 


1.48 


0.00 


S2 


224 


332 


1.48 


0.01 


S3 


228 


336 


1.47 


0.02 


S4 


229 


333 


1.45 


0.04 


S5 


230 


330 


1.43 


0.10 


S6 


230 


340 


1.48 


0.20 


S7 


233 


345 


1.48 


0.50 


S8 


235 


341 


1.45 


1.00 


Kl 


275 


418 


1.52 


0.01 


K2 


370 


560 


1.51 


0.01 


K3 


255 


404 


1.58 


0.01 



the various oscil l ation m odes. These modes were first pointed out 
bv IZanotti et aP diooi) in the investigation of the dynamics of 
purely hydrodynamical tori with nonconstant specific angular mo- 
mentum in Kerr spacetime, and are the consequence of the nonlin- 
ear coupling among modes, in particular of the / and oi modes. 

We complete our discussion of the spectral properties of these 
oscillating discs, by showing in Fig. |6] the PSD for model K2, 
which, we recall, represents a torus orbiting around a Kerr black 
hole with spin a = 0.7. As for the previous spectra, the dashed 
line corresponds to the unmagnetized version of model K2 and 
is included for reference. Overall, the features observed in a Kerr 
background are very similar to those found for models in the 
Schwarzschild case. Also in this case, in fact, the fundamental 
mode, its first overtones and the nonlinear harmonics are clearly 
identified and no evidence appears of new modes related to the 
presence of a toroidal magnetic field. 

As a final remark we note that the 2 : 3 ratio among the differ- 
ent p modes has a relevance also in a wider context. We recall, in 
fact, that among the several models proposed to explain the QPOs 



observed in LMXBs c ontaining a black hole candidate, the one sug- 
gested bv lRezzoUa et al. (2003a) is particularly simple and is based 
on the single assumption that the accretion disc around the black 
hole terminates with a sub-Keplerian part, i.e a torus of small size. 
A key point of this model is the evidence that in these objects the 
frequencies of the fundamental mode and the first overtone are in 
the 2 : 3 harmonic sequence in a very wide space of parameter. The 
simulations presented here further increase this space, extending it 
also to the case of magnetized tori and thus promoting the validity 
of this model for QPOs to a more general and realistic scenario. 



6.2 Gravitational-wave emission 

As pointed out bv lZanottietai] ( l2003h the oscillating behavior of 
perturbed accretion tori is responsible for significant changes of 
their mass quadrupole moment. As a result, these changes deter- 
mine the emission of potentially detectable gravitational radiation 
if the tori are compact and dense enough. This could be the case if 
the tori are produced via binary neutron star mergers or via gravita- 
tional collapse of the centr al core of ma s sive stars. In this section, 
we extend the analysis of IZanotti et al.l ( l2003l l2005h for unmag- 
netized discs and investigate the gravitational-wave emission from 
constant angular momentum magnetized tori orbiting around black 
holes. 

Although more sophisticated approaches involving perturba- 
tive techniques around black holes can be employed to study th e 
gravitational - wave e mission from these tori CNagar et al.l l l2005h : 
iFerrari et akl ( l2006h ; iNagar et all ( l2007h ). we here resort to the 
simpler and less expensive use o f the Newtonian quadrupole ap- 
proximation jZanotti et all ( l2003h ). which has been suitably mod- 
ifi ed to account for t he presence of a magnetic field, as done 
bv lKotakeetai] l l2004) . In particular, the quadrupole wave ampli- 
tude A20, and which is the second time derivative of the mass 
quadrupole moment, is computed through the "stress formula" 
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Figure 6. As Fig.|5]but for model ii'2. 



dObergaulinger et alj 1200^ ) 
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^20 



drdz 
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9$ r- 



(24) 

where k = 167r^/^/-\/T5, z = cos 6*, fij = pviVj — bibj, and $ is 
the gravitational potential, and is approximated at the second post- 
Newtonian order from the metric function grr = 1 — 2$ + 2$^. 

Figure [7] shows a spectral comparison between the designed 
strain sensitivity of the gravitational-wave detectors Virgo and 
LIGO, and the logarithm of the power spectrum \h{f)\y/J of the 
gravitational-wave signals for models S2, S8, and K2 (similar 
graphs are obtained also for the other models). Note that all of the 
sensitivity curves displayed in this figure assume an optimally in- 
cident wave in position and polarization (as obtained by setting the 
beam-pattern function of the detector to one), and that the sources 
are assumed to be located at a distance of 10 kpc. 

From Fig.|7]it is clear that all our models lie well above the 
sensitivity curves of the detectors for galactic sources and also that 
there are no significant differences in the power spectra as the mag- 
netic field strenght is increased. Interestingly, however, the signal 
from a torus orbiting around a Kerr black hole is clearly distin- 
guishable from the one around a Schwarzschild black hole. Besides 
having a fundamental mode at higher frequencies, in fact, also the 
amplitude is about one order of magnitude larger as a result of it 
being closer to the horizon and with a comparatively larger cen- 
tral density. As expected from the similarities in the dynamics, the 
signal-to-noise of these magnetized models is very si milar to one of 
the co rresponding unmagnetized tori, and we refer to lZanotti et al.l 
( l2005h for a detailed discussion. 



7 CONCLUSIONS 

We have presented and discussed the results of numerical simula- 
tions of the dynamics of magnetized relativistic axisymmetric tori 
orbiting in the background spacetime of either Schwarzschild or 



Figure?. Comparison between the power .spectrum v7 of the wave 

signal for models S2, S8, and K2 and the strain sensitivity of LIGO 
(dashed line), Virgo (long-dashed line). 



Kerr black holes. The tori, which satisfy a polytropic equation of 
state and have a constant distribution of the specific angular mo- 
mentum, have been built with a purely toroidal magnetic field com- 
ponent. The self-gravity of the discs has been neglected and, as the 
models considered are all marginally stable to accretion, the minute 
accretion of mass and angular momentum through the cusp is not 
sufficient to affect the background black hole metric. 

The use of equilibrium solutions for magnetized tori around 
black holes has allowed us to study their oscillation properties when 
these are excited through the introduction of small perturbations. In 
particular, by considering a representative sample of initial models 
with magnetic-field strengths that ranged from 2.5 x 10^"" G up to 
equipartition, and GRMHD evolutions over 100 orbital periods, we 
have studied the dynamics of these discs and how this is affected 
by a magnetic field. 

Overall, we have found the behaviour of the magnetized tori 
to be very sim i lar to t he on e shown by purely hydrodynamical tori 
jZanotti et al] ilOOj , |2005l) ). As in the hydrodynamical case, in 
fact, the introduction of perturbations triggers quasi-periodic oscil- 
lations lasting tens of orbital periods, with amplitudes that are mod- 
ified only slightly by the small loss of matter across the cusp. The 
absence of an inital poloidal magnetic field has prevented the devel- 
opment of the magneto-rotational instability, which could in fluence 
the os cillation properties and thus alter our conclusions ( jpragilel 
Determining whether this is actually the case will be the 
focus of a future work, where a more generic magnetic field con- 
figuration will be considered. 

As for unmagnetized tori, the spectral distribution of the 
eigenfrequencies shows the presence of a fundamental p-mode and 
of a series of overtones in a harmonic ratio 2:3:.... The analogy 
with purely hydrodynamical simulations extends also to the non- 
linear harmonics in the spectra and that are the consequence of the 
nonlinear coupling among modes, (in particular the / mode and 
of its first overtone oi). Also for them we have found a behaviour 
which is essentially identical to that found in unmagnetized discs. 
In summary, no new modes have been revealed by our simulations, 
and in particular no modes which can be associated uniquely to 
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the presence of a magnetic field. Nevertheless, the influence of the 
magnetic field is evident when considering the absolute values of 
the eigenfrequencies, which are shifted differentially to higher fre- 
quencies as the strength for the magnetic field is increased, with an 
overall relative change which is ~ 5% for a magnetic field near 
equipartition. 

Besides confirming the unmagnetized results, the persistence 
of the 2 : 3 ratio among the different p modes also has an im- 
portant consequence. It allows, in fact, to extend to a more gen- 
er al and realistic scenario the validity of the QPO model pres ented 
bv lRezzolla et alj j2003ah and lSchnittman &" Rezzollal |2003) , and 
which explains the QPOs observed in the x-ray luminosity of 
LMXBs containing a black hole candidate with the quasi-periodic 
oscillations of small tori near the black hole. The evidence that this 
harmonic ratio is preserved even in the presence of toroidal mag- 
netic fields, provides the model with additional robustness. 

When sufficiently massive and compact, the oscillations of 
these tori are responsible for an intense emission of gravitational 
waves and using the Newtonian quadrupole formula, conveniently 
modified to account for the magnetic terms in the stress-energy ten- 
sor, we have computed the gravitational radiation associated with 
the oscillatory behaviour. Overall, we have found that for Galactic 
sources these systems could be detected as they lie well within the 
sensitivity curves of ground-based gravitational-wave interferome- 
ters. 

As a concluding remark we note that our discussion here has 
been limited to tori with magnetic fields whose pressure is at most 
comparable with the gas pressure, i.e., j3c < 1. The reason behind 
this choice is that while in equilibrium, the magnetized tori with a 
purely toroidal magnetic field are not necessarily stable. Rather, in- 
dications coming both perturbative calculations and from nonlinear 
simulations, suggest that these tori could be dynamically unstable 
for sufficiently strong magnetic fields. The results of these inves- 
tigatio ns will be presented in a forthcoming paper jRezzoUa et al.l 
1^03)). 
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